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ABSTRACT 



In a recent paper, ISalvador-Sole et al.1 (|2012f) have derived the typical inner structure 
of dark matter haloes from that of peaks in the initial random Gaussian density field, 
determined by the power-spectrum of density perturbations characterising the hierar- 
chical cosmology under consideration. In the present paper, we extend this formalism 
to the typical kinematics and triaxial shape of haloes. Specifically, we establish the 
link between such halo properties and the power-spectrum of density perturbations 
through the typical shape of peaks. The trends of the predicted typical halo shape, 
pseudo phase-space density and anisotropy profiles are in good agreement with the 
results of numerical simulations. Our model sheds light on the origin of the power- 
law-like pseudo phase-space density profile for virialised haloes. 
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1 INTRODUCTION 

Virialised haloes in iV-body simulations of cold dark matter (CDM) cosmologies show a wide variety of ellipsoidal shapes. 
On the contrary, their structural and kinematic properties are remarkably similar from one object to another. They are little 
sensitive not only to the mass, redshift, environment and even specific cosmology, but also to their individual shape. Only 
their scaling shows a mild dependence on some of these properties. As shown by gravo-hydrodynanic simulations, baryons 
introduce a larger scatter in the properties of haloes at their central region. However, in the present paper, we will concentrate 
on pure dark matter haloes and we will not deal with such secondary baryonic effects. 

The typical spherically a veraged halo density profile, (p)(r), is well-fitt e d, dow n to about one hundredth the virial radius, 
by the so-called NFW profile jNavarro et al.ll 1997^) as we l l as by the Einasto (1 19651) profile, which give s slightly better fits down 



to smaller radii (jNavarro et al]|2004l ; iMerritt et al]|2005l . I2OO6I : Istadel et"al]|2009l ; [Navarro et alboiol ). The velocity disp 



profile, cr(r), is reason ably well-fitted by the solution of the J eans equation for sp herically symmetri c isotr opic systems with 
null value at infinity I Cole fc Lacev 19961 : Merritt et al. 20061 ). More remarkably, Taylor fc Navarrol ( 2001 ) showed that the 
pseudo phase-space density profile is very nearly a pure power-law, 

(P)(r 



= Ar v . 



(1) 



with index v : 



1.875 (see also Ascasibar et al. 2004; Rasia et al. 2004: Bctt ct al. 2007; Faltcnbacher et al 



2007 



Vass et al 



2009; [Navarro et afcoiOT ) and a similar relation arises from the radial velocity dispersion component, a T (r) (jDehnen fc McLaughlin! 
2005). In equation {TJ, a(r) is the velocity dispersion over the whole spherical she ll with r, so the varianc e coincides with 
the spherical average of the local value at points with r, <r 2 (r) = {o~i oc )(r). Finally, Hansen fc Moorel ( 20061 ) found that the 
velocity anisotropy profile /3(r) behaves linearly with the logarithmic derivative of the density, 



p(r) = 



dln(p) 
dlnr 



+ b 



(2) 
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with a and b respectively equal to about —0.2 and 0.8 ( Hansen fc Stadel 20061 ; see also Ludlow et al. 2011 for an alternative 
expression), although with a substantial scatter this time. 

The origin of all these trends is certainly related with the way dark matter clusters. In hierarchical cosmologies, haloes 
grow through continuous mergers with notably different dynamic effects according to the relative mass of the captured 
and capturing objects. For this reason, it is usually distinguished between major mergers, with a dramatic effect each, and 
minor mergers, contributing together with the capture of diffuse matter (if any) to the so-called accretion, responsible of 
just a smooth secular evolution of the system. Some aut hors have attempted to explain the typical halo density profile as 



the result of repeated m ajor (or intermediate) mergers ( Sver fc White! E.998I ; ISalvador-Sole et al. 19981: 



Subramanian et al. 



2000l; iDekel et al.ll2003l). Others have concen t rated instead in the e ffects of pure accretion (PA) IIAv ila- Rees e et al 



1998; 



Nusser fc Shethl Il999l ; bel Popolo et al.l I2OO0I ; llVIanrique et all I2OO3I ; Uscasibar et all |2004| ; ISalvador-Sole et al.ll2007f h Both 



extreme scenarios have also been investiga ted regarding the possible origin of the pseudo phase-space dens ity and velocity 
anisot ropy profiles ( Hansen fc Moore 20061 ) . The PA scenario has received much support from the results by Wang fc White! 
(2009) showing that all typical halo trends are already set in the first generation haloes formed by monolithic collapse (i.e. 
no major merger; only accretion of diffuse matter) in warm dark matter cosmologies. 



(e-K 
200( 



Regarding the shape, CDM haloes are found to be triaxial ellipsoids, with a trend towards prolate rather than oblate shapes 
iFrenk et allll988l; IPubinski fc Carlb"er3ll99ll; Iwarren et al.ll 19921; ICole fc Lacevlll996l; ISpringel et al]|2004l; Ullgood et Z 



2006; lHavashi et al.l 120071 ; IMaccio et al.l |2007; Stad el et alj 120091 ; IVera-Ciro et al.1 



2011 



). Inside each individual object, the 

typical mi nor to major axial ratio takes a roughly uniform value of about 0.6, w i th a slight trend to an outward-decreasing 



triax i alitv jFrenk et allll988l;lBullockll2002l;|jing fc Sutoll2q02l; ISpringel et al.ll20o4 iKasun fc Evrardl 



20051; Ullgood et alj|2006l ; lHavashi et al.ll2007l ; iBett et al.l 120071 : IStadel et al.l I2OO9I; IVera-Ciro et al 



20051; Bailin fc Steinmet3 



2011). The main axis is 



Lemson fc Kauffmannl 19991 ; Basilakos et al 



preferentially aligned, at all scales, along with the filament feeding the halo (e.g. ^_ 

20061; IPatiri et allbooej ; IMaccio et al1l2007l ; iRagone-Figueroa et allboiol ; IVera-Ciro et al.ll201ll ). Th is indicates that the 



ory of the preferred direction of major mergers and accretion is not erased during virialisation (|Vera-Ciro et al.l l2011f ) or, 
equivalently, that the shape of virialised haloes depends on that of their seeds. Moreover, as haloes are not supported by 
rotation but by the local anisotropic velocity tensor, the fact that their triaxial shape is related to the shape of protohaloes 
automatically implies that their kinematics must also be related to it. 

The seeds of haloes are believed to be peaks (secondary maxima) in the primordial rando m Gaussian density field filtered 
at the scale of the halo. The isodensity contours in the immediate vici nity of peaks are triax ial i|Doroshkevichlll970h and rather 
prolate with a trend to become more spherical for very high peaks dBardeen et al. 19861 ; hereafter BBKS). As well-known, 
the monolithic collapse of non-spherical systems is highly non- radial (|Zeldovicb]|l97oh~ giving rise to filaments and triaxial 
virialised objects. Thus, it is natural to believe that th e shape of peaks is somehow transl ated into that of haloes, in agreement 
with the above mentioned alignments. A few authors ( Lee et al. 20051 : Rossi et al.ll201 ll ) have tried to make the link between 
the shape of haloes and that of peaks through the modelling of ellipsoidal collapse. Unfortunately, these models do not account 
for the highly non-linear effects of shell-crossing during virialisation, which play a crucial role in setting the final properties 
of virialised haloes. On the other hand, there is in the literature no attempt to relate the kinematics of haloes with the shape 
of their seeds. 

In a recent paper, ISalvador-Sole et al. (2012: hereafter SVMS) have shown that the kinematics of virialised haloes in 
(bottom-up) hierarchical cosmologies with dissipationless collisionless dark matter depends on their triaxial shape, contrarily 
to their spherically averaged density profile which does not. This allowed SVMS to infer, under the assumption of PA, the 
typical spherically averaged density profile for haloes from that of peaks in the primordial density field, determined by the 
power-spectrum of density perturbations. Furthermore, SVMS showed that the density profile for haloes having undergone 
major mergers is indistinguishable from that of haloes grown by PA, so the model actually holds for all haloes regardless of 
their individual aggregation history. 

In the present paper, we extend the SVMS model to the kinematics and triaxial shape of virialised objects. Under 
the PA assumption and neglecting any possible rotation tidally induced by surrounding matter, we derive the halo shape, 
velocity anisotropy profile and pseudo phase-space density profile from the triaxial shape of peaks, taking into account the 
full virialisation process. We first assume the simple case of PA and then study the foreseeable effects of major mergers. This 
allows us to establish the link between those typical halo properties and the power-spectrum of density perturbations. The 
theoretical predictions obtained when this formalism is applied to CDM haloes are in good agreement with the results of 
numerical simulations. 

The paper is organised as follows. In Section [2] we derive some general relations valid for triaxial systems, regardless of 
whether they are in equilibrium or not. Assuming PA, these relations are used, in Section [3] to make the link between the 
triaxial shape of a virialised object formed by PA and that of its seed. The typical velocity anisotropy and velocity dispersion 
profiles for virialised objects are derived in Section [4] from their triaxial shape. In Section [5] we apply the model to CDM 
haloes. The origin of the power-law-like form of the pseudo phase-space density profile and the effects of major mergers are 
discussed in Section [6] Our results are summarised in Section [7] A package with the numerical codes used in the present paper 
is publicly available from www. am.ub. es/~cosmo/haloes&peaks . tgz. 
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2 MEAN SQUARED AND CROSSED FLUCTUATION PROFILES 

There are in the literature several ways to characterise the shape of a triaxial sys tem with semiaxe s a ^ b c: the axial 
ratios, the ellipticity and prolateness (e.g. BBKS) and the triaxiality parameter fe.g. iFranx et aUll99lT ). among others. In the 
present paper, we use the primary and secondary eccentricities, respectively defined as 



e P = 1 ~ — 

or 



1/2 



and 



or 



1/2 



(3) 



As the eccentricities vary, in general, over the system (we assume from now on that the centre of symmetry remains the same 
at all scales), we will deal with the eccentricity profiles c p (r) and e s (r), being r the radius of a sphere associated with each 
ellipsoid (see details below). 

Numerical studies of virialised haloes usually work with spherically averaged profiles. Of course, triaxial systems show 
non-null departures 5p(r,9,(p) and <5$(r, 6,<fi) of the local density p(r,9,ip) and gravitational potential <6(r, 6, if) from their 
respective spherical averages (p)(r) and (5>)(r). Hence, the (scaled) mean squared or crossed density and potential fluctuation 
profiles, {(Sp/ (p)) 2 ){r), ((5$/(<3?)) 2 )(r) and ((5p/(p))(5$/(<3?)))(r-), also give a measure of the triaxial shape of those systems. 
Below, we relate these profiles to the eccentricity profiles defined above. 

Take the Cartesian axis z aligned along with the major axis of the ellipsoidal isodensity contour p lso (r), labelled by the 
radius 

1/2 



5 (a + b 2 + , 



The local density at the point (r, 8, <p) then takes the form 



p(r,6,tp) = p 1BO (r) 



Sp(r) + ef(r) 



sin 2 8 cos 2 c 



■ + 



+ 



efp(r) 1 - ef s (r) 



(4) 



(5) 



where e ps stands for one of the two eccentricities and e sp for the other onqj. The spherically averaged density at r is then 
equal to 

eg(r) + e s » 



(p)(r) = 



P ls °(r) 



1 



where we have defined the function 
1 1 



G(r) = 



1 + 



l-e 2 (r) 



+ 



4(r) 



G(r) 



= a (r) 



1 1 1 

+ TT7 — r + 



o 2 (r) b 2 (r) c 2 (r) 



Dividing equation © by {p)(r) and replacing p lso (r) by its value given in equation we obtain 



Sp , 

1 + V) ir ' 



sin 2 8 cos 2 6 + 



+ 



G(r) [ ~ ^ ' l-e 2 p (r) ' l-e 2 s (r)_ 
The mean squared density fluctuation over the sphere with radius r, 

2\ rir (.Itx / „ \ 2 



dO sin 8 



( r ,e,< 



can then be readily obtained from Sp/(p) given by equation ([8|. The result is 




5 1 [(f-e 2 )(f- e 2 ) + (l-e 2 ) + (l-e 2 )] 2 



:) 2 + (l^e 2 ) 2 + (l-e 2 ) 2 ] 



(r). 



(6) 



(7) 



(8) 



(9) 



(10) 



The relations between the eccentricity profiles and the mean squared potential and crossed density-potential fluctuation 
profiles follow from the relation (|10p and the relations between such fluctuation profiles and the mean squared density 
fluctuation profile. To derive them we need first to develop the member on the left of the Poisson equation, 



(*) 1 + 



(5$ 

W> 



= 4ttG(p) 1 + 



5p_ 

(P) 



(11) 



Taking into account that, by the Gauss theorem, ($) satisfies the usual Poisson integral relation for spherically symmetric 
systems 



1 The specific values of e ps and e sp depend on the orientation of the x and y Cartesian axes relative to the minor and intermediate 



semiaxes. 
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d($)(r) _ GM(r) 

dr r 2 ' 

we arrive after some algebra at the exact relation 



Sp_ s$_ _ M d_ ( 5$ 

(p) {$> ^ r)T ' dr \ <$> 



^(r)C(r)r 2 V 2 



5$ 



(12) 



(13) 



where £(r) = (dlnM/dlnr) -1 and £(r) = — (din | ($) |/d In r) _1 . Equation (fT3)l shows that it is possible to infer 5$/(<&) 
from 5p/(p) for the appropriate boundary conditions. This requires, however, the full characterisation of 8p(r,0,(p), which 
is, in general, a complicate function of the spatial coordinates. For simplicity, we will concentrate in ellipsoidal systems 
(i) with non-rotating symmetry axes and (ii) with rms density fluctuation profile of the power-law form, ((Sp/ (p)) 2 } 1 ^ 2 = 

nr~ 1 5p(r, 8,tp)/ (p)(r)|l Condition (i) is essentia l ly satisfied by 



simulated haloes (see 


Jhm & Suto 2002; Springel et al. 2004 


2007: Vera-Ciro et al. 


2011), while condition (ii) can be seen 


rang^E 



In these conditions, multiplying equation (|13[1 by the scaled density fluctuation, performing the spherical average and 
taking into account the divergence theorem, we are led to 




K{r) 



2 [1 + (1 - «)C(r)] d 



CM 



dr 



Sp (5$ 

Tp)W> 



+ 



dr 2 



Sp 54> 

MM 



(14) 



being K(r) = 1 + K^(r)[2 + (1 — k)£(t)]. Equation (|14[) is an ordinary differential equation that can be solved, for trivial 
consistency boundary condition at r = 0, for the mean crossed density-potential fluctuation. Similarly, multiplying equation 
()13[) by the scaled potential fluctuation, the same derivation as above leads to 



M (*) / U 




CMC(r)r 2 




(15) 



Equation (|15p cannot yet be solved for the mean squared potential fluctuation because of the presence of the unknown mean 
squared radial derivative of the scaled potential (the last term on the right). But this drawback can be bypassed by multiplying 
equation (1 13 p by the partial derivative of the scaled potential fluctuation and operating as usual. This leads to 



d_ / 5p_ 

dr\(p) <$) 



i^V) 




<(r)({r)r ■ 



2[1 + C(r) 
CM 



d S$ 
dr ($) 




d 8JT 
dr'J®). 



(16) 



Then, substituting in equation (|16[) the mean squared partial derivative of the scaled potential given by equation (|15|l and its 
radial derivative drawn from the differentiation of the same equation, we arrive at 



Sp <5<E> \ 3 r d 

M W/ (r)+ 4iV)d^ 



5p_ te\ = £(r) S± ] )lr) 

(p)<*>/ 7m\VW. ; 



•Hfet J(r)A 
4J(r) I dr 




+ [1 + CM]^ 




being 

JM = 



4/(r) + 3CM-2r 1 M + ^ + 5 
dmr 



J(r) = 1 + C _1 (r) + 



1 din 
4 dlnr 



_C_ 

?r 2 



(17) 



(18) 



and 7(r) = J(r) — «;/2. Once again, equation (|17[) can be solved, for trivial consistency boundary conditions at r = 0, for 
the mean squared potential fluctuation. Note that the mean d(Sp/ {p}) /dr <5<J>/{$) profile entering the Jeans equation for 
anisotropic spherically averaged triaxial systems (eq. [42] below) can also be inferred from the radial derivative of the mean 
crossed density-potential fluctuation and the condition d[Sp(r,9, tp)/ (p)(r)]/dr — Kr~ 1 Sp(r,6, <p) / (p)(r) above. 

One particular case of interest is that of homologous triaxial systems, i.e. with ft = 0. Equations (|14[) . (|15|l and (|17[) then 
have the trivial solution, 



2 Such an equivalence can be proven taking into account that, as {Sp/ (p)) = 0, either the integral of 5p(r, 9, <p)/ (p)(r) over cj) is an even 
function of 9 or the integral of svn(9)Sp(r, 9, ip)/(p) (r) over 9 is an odd function of <f> and that the condition ({Sp/ (p)) 2 } 1 / 2 = Qr K implies 
the contrary relation for {d[Sp(r, 9, ip)/ {p)(r)]/dr} — n[Sp(r, 9, ip) / (p)(r)]/r. 

3 In fact, an arbitrarily accurate approximation can be achieved by splitting the real log-log mean square density fluctuation in a series 
of concatenate linear functions so that such a condition would hold in each small segment. 
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Figure 1. Predicted Q/P ratio (see the colour bar for the different ranges of this quantity) between the rms density and potential 
fluctuation profiles in self-similar objects with k the power index of the rms density and potential fluctuations and —a that of the 
spherically averaged density profile. 



Sp (5$ 



(r) = 




(r) = 




(r) = Q. 



(19) 



That is, all squared or crossed fluctuation profiles are uniform and equal to each other. 

Another particular case of interest is that found in asymptotic regimes, i.e. with all the profiles behaving as power-laws. 
Equations (|15|) . (|16|) and ()17p then imply that both the rms density and potential fluctuation profiles have the same form, 
proportional to r K . This does not mean that the isodensity and isopotential contours are equally non-spherical at any given 
radius: the positive proportionality factors Q and P of the rms density and potential fluctuation profiles, respectively, are 
different in general. Specifically, for an object with power-law spherically averaged density profile (p)(r) oc r~ a , the ratio Q/P 
depends on k, and a in the following way, 



1 - 



(1 — o)(k — 1)k 
(2-q)(3-q) 



1 - 



("-!) 
(2-a) 



1 

K— 1 



6 k (3 — a) 



2a + k/2 



(20) 



as inferred from equations ()14[) and (|17|) . Figure [T] shows the ratio Q/P that results for a and k respectively in the ranges 
2 < a < 3 and —2 < k < 2. The function is well-behaved everywhere except in the line k = 4a — 10 where it diverges by 
alternating positive and nonphysical negative values on opposite sides of the line, in a series of segments separated by a few 
well-behaved points. Far enough from the strip, Q/P tends to take values greater than one for any value of k. 

Notice that, as p(r,0,<p) is positive, the negative values of Sp necessarily satisfy the condition \Sp\ < (p). Moreover, in 
rather prolate ellipsoidal systems the solid angle where Sp is negative is less than the solid angle where it is positive and, as the 
mean Sp is null, we also necessarily have that the positive values of Sp are less than (p). Consequently, we have \5p/{p)\ < 1, 
in general. Then, provided the mass distribution is smooth enough, which is guaranteed in triaxial virialised objects, equation 
(|13[) also implies < 1. The angular cross-correlation of these two dimensionless quantities is even smaller, so any 

physical quantity reporting to the triaxial system can be written as a series expansion of the dimensionless fluctuation profiles 
measuring the deviation from spherical symmetry. We will take advantage of this fact in Sections [3] and [4] 

All the previous relations hold regardless of whether the triaxial system is in equilibrium or not. Thus, they apply both 
to the virialised object and its seed. From now on, all quantities referring to the protoobject are denoted by subindex p. 



3 ECCENTRICITY PROFILES 

We want to relate the eccentricities of a virialised triaxial object with those of its seed. As there are two eccentricities, we 
need two equations. These equations must obviously report on properties involving the shape of these systems. 
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One of such properties is the volume of ellipsoidal isodensity contours. The ratio between the volumes of the ellipsoids 
encompassing a given mass in the final and initial systems must be equal, to leading order in the deviation from spherical 
symmetry, to the ratio between the volumes of the corresponding spheres encompassing identical mass, 

(21) 



a(r) b(r) c(r) 



a p (rp)6p(rp)cp(r p ) rj(r) ' 

To write the member on the right of equation (|2ip we have taken into account that the mass of an ellipsoidal isodensity 
contour with semiaxes, a, b and c, labelled by R 

M{R)= I AO sine I d(j> drr 2 p(r,6,(f>), (22) 
Jo Jo Jo 

coincides, owing to the relations ([5j and |6| arising from the definition Q of the labelling radius, with the mass of the sphere 
with radius R, 

r-R 



M(R) = 4tt / drr (p)(r). 
Jo 

On the other hand, the function r p (r) in equation pip is the solution, for the boundary condition r p (0) 
equation 

dTp = r 2 (p)(r) 
dr r% (p P )[r p (r)} ' 



(23) 



0, of the differential 



(24) 

that follows from differentiation of equation (|23[) holding for spheres with identical mass in the initial and final systems. For 
simplicity in the notation, we omit from now on the explicit dependence of r p on r. 

We stress that equation (|2ip is only valid to leading order in the deviations from spherical symmetry. This does not mean 
that the volume of each ellipsoid is approximated by that of the corresponding sphere. (By doing this, we would lose the 
information on the shape of the system.) What we approximate is the whole ratio between the volumes of both ellipsoids by 
the ratio between the volumes of the corresponding spheres. Were the axial ratios of the ellipsoid conserved over the evolution 
of the system, the relation (|21[) would be exact. Actually, the axial ratios vary during virialisation. But this variation is of 
higher order in the deviation from spherical symmetry, so the ratio between the volumes of the two ellipsoids is kept equal to 
that between the two spheres to leading order. 

Taking into account the relation Q, equation (|21jl leads to 



4)0- 



(1-egXl-eg) 
.(l-eg) + (l-e 



(25) 



[l + (l- e |) + (l-el)]3w + 

A second property involving the triaxial shape of the system is the energy of the sphere with radius R encompassing a 
fixed mass M (see SVMS), 



E{R) = Aty drr 2 {p)(r) 



a 2 (r) GM(r) 



2tt 



dr r 2 (SpS&)(r) 



This can be written in the form E(R) 
£{R) = 4x [ drr 2 (p}{r^ 



S(R) + 5£(R). The so-called "spherical" total energy, 
GM{r) 



(26) 



(27) 



measures the total energy in the sphere with M, were the mass distribution inside it spherically symmetric and endowed with 
the velocity variance profile s 2 (r) that would result taking into account the energy lost by shell-crossing during virialisation 
but not the potential energy exchanged with the rest of the system owing to its non-spherical mass distribution. In the 
protoobject, Sp(r p ) can be taken equal to the real velocity variance profile <Tp(r), but in the virialised object s 2 (r) differs from 
u 2 (r). Nonetheless, the profile s 2 (r) can be obtained, by differentiation of equation (1271) . 



d£/dr + GM(r) 



S{r) = 2 dM/dr 
from the mass and spherical energy profiles, M(r) and £(r) 



respectively, derived in SVMS. The result was 



r{M) = 
and 

£(R) = 



3 GM 2 



10 £ P (M) 



-R 



dr 



GM(r) dM + W(r) 



dr 



2r 2 



(28) 



(29) 



(30) 
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where we have defined the "spherical" potential energy, 



W(R) = -4tt 



drr 2 <p>(r)^ 



This leads to the following explicit expression 



s 2 (r) 



47rr 2 (p)(r) 



W(r) 



df 



87r(p)(f)GM(f) + 



W(f 



(31) 



(32) 



showing that s 2 (r) can be calculated from the spherically averaged density profile (p)(r). As this latter profile is independent 
of the shape of the system, so must the function £ (r) (see eq. [27]). Consequently, all the effects of triaxiality in E(r) are 
included in the residual 



8£ (R) = 2tt / dr r (p) (r) ($) (r) 



5p 5$ \ 
<P> <<&> / 



(r) + 



a 2 {r)-s 2 (r) 



<*)(r) 



(33) 



Specifically, the part of it harbouring the difference a 2 (r) — s 2 (r) corrects the spherical total energy in the sphere from the 
gravitational energy exchanged (it may be positive or negative) at a distance with the rest of the system, and the other 
part with the crossed density-potential fluctuation profile corrects it for the different potential energy due to the actually 
non-spherically symmetric mass distribution of the system. 

As virialisation is driven by the energy lost through shell-crossing (not by the energy exchanged at a distance between 
shells, which is a consequence of the former effect and much less important) , the ratio between the final and initial total energies 
in the sphere is, to leading order in the deviation from spherical symmetry, equal to the ratio between the corresponding 
spherical total energies also accounting for that energy loss 

E(R) £{R) 



E P (R P ) £ P (R P ) 



(34) 



We stress, once again, that equation (|34|) is valid to leading order in the deviation from spherical symmetry. This does 
not mean that the total energy in the sphere is approximated by its spherical counterpart. (By doing this we would lose all the 
information on the shape of the system.) ft is the whole ratio between the total energies in the final and initial spheres which is 
approximated by the ratio between their spherical counterparts. Thus, taking into account the relation E(R) — £(R) + 5£(R), 
equation (|34|) implies 

(35) 



£(R) _ S£(R) 
£ P {R) ~ S£ P (R) ' 

to leading order in the deviation from spherical symmetry. 

Substituting the expressions for 8£ and 8£ p in the virialised object and its seed (eq. [33] 
differentiating it, we are led to 

a 2 (r)-s 2 (r) ] I 8p p 8$ p 

wr =<$p>(rp) W(^y, 



into the relation (1351) and 



WW 

VM 



8p 8® \ 

V)W) 



(r) + 



where T>{R) = £(R)/£ P (R), the so-called "spherical" energy dissipation factor, is given by (eqs. [29] and [30| ) 



V(R) 



WR 2 



3GM 2 (R) 



dr 



47r(p)(r)GM(r) + 



W(r) 
2r 2 



(36) 



(37) 



Then, taking into account equations (|10[) and (|14[l . both for the virialised object and its seed, equation (|36[) establishing the 
link between the respective mean crossed density-potential fluctuation profiles leads to 

3[(l- e 2 ) 2 (l- e 2 ) 2 + (l-e 2 ) 2 + (l-e 2 ) 2 ] 



U(r 
where 
U(r) = 



3[(1 



4) 2 (1 



+ (!• 



+ (1 



[(l-e 2 )(l-e 2 ) + (l- e 2 ) + (f 



SUr)-. 



[(l-e 2 )(f - e 2 ) + (l-e 2 ) + (l-e 2 )] 2 



(* P )(r p )2>(r)V(r) 



and 



Sir) 



5a 2 (r) 
2 



('•) 



(*>(r) 



V(r). 



(r p ), (38) 



(39) 



being V(r) — 1 — C( r )7( r ) {l — [1 + 2t]7(r) — d]n ^ | and 7(7-) the logarithmic derivative of the mean crossed density-potential 
fluctuation profile. 

Equations (|25[) and ()38[) determine the eccentricities of the virialised object grown by PA from those of its seed, as 
wanted. From equations (|38[) and (|39|) . we see that this relation involves the profile <r 2 (r), which depends on the shape of the 
final system and, hence, on that of its seed. But such a dependence cannot be determined without previously determining 



© 0000 RAS, MNRAS 000, 000-000 



8 Salvador-Sole, Serra, Manrique & Gonzdlez-Casado 



i 




0.8 



0.6 



a 



<D 



0.4 



0.8 











1 



2 



3 



4 



U(R) 



Figure 2. Primary eccentricity space (delimited by the solid red line) possible to occupy at small and intermediate radii by current 
CDM haloes with arbitrary mass in the concordance model for varying values of the quantify U(R) (see text). The two horizontal black 
lines delimit the values of e p reported in the literature for simulated haloes. 

the anisotropy profile. This will be done in Section [4] Only at small r, the term S(r) in equation (|38|) becomes negligibljf]. 
This means that, except for the factor U(r), equation (|38[) takes there the form of an identity relation, just as equation (|25l) . 
Interestingly, the set of algebraic equations (|38[) and (|25[) is solvable only for a very narrow range of U(r) values around unity. 
This is shown in Figure[2l where we plot the solution space for varying values of U(r) inferred by means of the algebraic solving 
procedure provided in the Appendix. Consequently, the eccentricities of the virialised object near the centre are necessarily 
close to those of the protoobject at the same location. 

4 KINEMATIC PROFILES 

As mentioned, the kinematics of a non-rotating triaxial virialised object must depend on its shape and conversely. This 
relation is however hard to establish from first principles, even under the assumption of PA as made here. It can nonetheless 
be guessed from comparison with the simplest idealised case of strict spherical symmetry. See also Section [5] for a more 
fundamental justification. 

In that idealised case, the final velocity dispersion would be purely radial because the system would expand, collapse and 
virialise radially (shell-crossing would take place between spherical shells). The total energy within spheres is, to leading order 
in the dimensionless fluctuation profiles measuring the deviation from spherical symmetry, the same in spherically symmetric 
systems as in triaxial ones (see eq. [26]). Thus, non-radial infall due to non-spherical shells (or to any possible instability; 
see below) should cause the transfer from radial to tangential kinetic energy together with a deviation of the gravitational 
potential from its spherical average without altering, to leading order, the total energy of the system. Assuming one-to-one 
correspondence between the two effects, we are led to the conclusion that, to leading order, the fractional velocity variance 
transferred from the radial to the tangential direction (equal to half the fractional 1-D tangential velocity variance generated) 
should be equal to half the typical fractional deviation of the potential from its spherical average, 



Equation (|4UI) might seem to imply that it is impossible to have tangential velocities in spherically symmetric virialised 
objects (with null (5$), which would be manifestly wrong. What it actually implies is this very conclusion but for objects 

4 The cr 2 (r) profile at every radius is of the order of the squared circular velocity, GM(r)/r, while becomes much larger than 

GM(r)/r at small radii. 




(40) 
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formed by Indeed, as discussed in SVMS, such objects are always ellipsoidal (they have non-null 5$). Even if the seeds 
are spherically symmetric, the gravitational effect of the surrounding matter (non-null shear tensor) or radial orbit instability, 
in the artificial case of isolated seeds, automatically leads to non-radial collapse and the formation of triaxial virialised objects 
(see Sec. [5] for more details). 

Given the relation (|40|l . the anisotropy profile, 



/3(r) 



(r) 



?(r) 



1-2 



(41) 



is a function of the rms potential fluctuation profile for the virialised object. As shown in Section [3] such a fluctuation profile 
is related to its counterpart in the seed. However, such a relation involves the velocity dispersion profile itself (eq. [38]). Thus, 
tj 2 (r) must be determined before the rms potential fluctuation profile. 

To do this we will make use of the (exact) generalised Jeans equation for steady non-spherical virialised objects (see 
SVMS), which, taking into account equation (|41|) . can be written as 



d[(3-2/?)- 1 ( /0 >a 2 ] , 2/3(r) (p)(r)a 2 (r) 



dr 



• + 



2p{r) 



+ (p)(r) 



GM{r) 



= -<p)(r)(#>(r) 



2 / 5p <5$ \ 



r \ (p) ($) / 



(p) ) <<£-> / 



(42) 



The member on the right of equation (|42p is of second order in the dimensionless fluctuation profiles. Like fl(r), these profiles 
are functions of their counterparts in the seed and of o~(r). Consequently, equation (|42[) is a differential equation for <r(r), 
which can be solved for the usual boundary condition of null dispersion at infinitjjf]. Moreover, the member on the right of 
equation ()42|) is of second order in the deviations from spherical symmetry, while all remaining terms are of lower order; the 
anisotropy profile itself is correct only to first order. Therefore, it can be neglected in equation (|42p . which notably simplifies 
the solution of this equation. Once <r(r) has been determined, the eccentricities of the virialised object can be calculated from 
equations (|25[) and (|38|l and /3(r) from equation (|41[) . 



5 APPLICATION TO CDM HALOES 

In SVMS, the spherically averaged density profile for CDM haloes was shown to arise from the spherically averaged density 
profile of halo seeds. As this latter profile is determined by the power-spectrum of density perturbations, it was possible to 
infer the typical halo density profile for CDM haloes directly from the power-spectrum of density perturbations. The extension 
of the model presented in the present paper relates the shape and kinematics of CDM haloes to the shape of their seeds, 
which is also determined, of course, by the power-spectrum of density perturbations. Thus, it should also be possible to infer 
the typical shape and kinematics of haloes directly from the power-spectrum of density perturbations. 

To do this we need to know the eccentricity profiles for protohaloes, while what is only known is the typical eccentricities 
at peaks (i.e. at r p = 0) in the initial density field filtered by a Gaussian window, calculated by BBKS. Nonetheless, the former 
can be accurately inferred from the latter, following the procedure developed in SVMS for the spherically averaged density 
profile. 

As pointed out in SVMS, the fact that haloes grown by PA develop from the inside out implies that the density profile 
for each halo ancestor in the continuous series leading to the final object exactly matches that for the halo within one triaxial 
isodensity contour. All these halo ancestors evolve from peaks in the initial density field filtered at the ancestor scale. Thus, 
the typical spherically averaged density contrast profile for the protohalo, (S p )(r p ), convolved by a Gaussian window with any 
radius Rf must be equal to the typical density contrast of a peak, S p ^(Rf), in the density field filtered at scale Rf, 



1/2 



6 Bk (Rt)={~) i? f - 3 J dr p r 2 (5 p )(r p )exp(^^|j . (43) 

The typical trajectory <5 p k(-Rf) followed by pea ks evolving by PA into the series of ancestors leading to a typical halo with M 
at t is the solution of the differential equation ( Manrique fc Salvador-Sole! 1996T ). 



^ = -Xe(Rf,6 pk )0-2(Rl)Rt , (44) 

for the appropriate boundary condition (through the <5 p k(£) and Rf(M) relations given by eqs. [40] and [41] of SVMS). In 
equation (|44p. the inverse of x c (Rf, S pk ) is the mean inverse curvature x (i.e. minus the Laplacian over the second order spectral 



5 In Section [6] we will see that it also holds for virialised objects having undergone major mergers, although, in this particular case, 
such an implication is quite obvious. 

6 The inside-out growth of object formed by PA guarantees that the solution out to any radius is kept unaltered as the radius of the 
object increases, so we have indeed the right to endorse the boundary condition at infinity. 
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moment) of peaks with <5 p k at Rf. The differential equation (|44|) can be solved, which then allows one to inverse the Fredholm 
integral equation (|43[) and obtain the wanted density contrast profile for the typical protohalo (see SVMS for details). This 
procedure was followed in SVMS to infer the unconvolved spherically averaged density profile for the seed of a typical halo. 

In the present case, we must take into account that the squared semiaxes of each peak over the peak trajectory leading 
to the typical halo with M at t coincide (except for an arbitrary scaling) with minus the second order spatial derivatives \? 
of the filtered density contrast at the peak, in Cartesian coordinates Xj with xi aligned along the major semiaxis. They are 
therefore related to the density contrast profile, 5 p (r), for the protohalo through 



(45) 



Xf( Rl ,5 pk )a 2 ( Ri ) = jdr p( 5 p (r p ) exp [-± p^)"] J^- 

Writing the density contrast profile in the halo in terms of its spherical average (this relation is similar to that between the 
respective total densities, given by eqs. [5] and [6]), performing the second order spatial derivatives of the Gaussian window 
inside the integral on the right of equation (|45[1 and taking x = 0, changing to spherical coordinates, integrating over and 
<f) and averaging over the peak curvature, we arrive at 



If * 



U - ' j I {Rf, <5 P k) <T 2 (i? f ) R 2 { 



J drprp ^ p ^ rp ^ exp \ 2R 2 { 



r? 



SH?y drp ^ (5p>(rp) ^ j( S P p V P ) eXP l"Wj' (46) 
where {{x))(R{, <5 p k) is the mean curvature of peaks with <5 p k and Rf and A,(Rf, 5 p k) (j = l-r-3) are its associated dimensionless 
semiaxes Aj scaled to the square root of the Laplacian, so that they satisfy A\ + A% + A\ = 1. These semiaxes depend on <5 p k 
and f?f through {{x)) in a well-known form calculated by BBKS. Thus, we can calculate them as well as the curvature {{x)) 
over the peak trajectory <5 p k(-Rf) leading to the halo with M at t and invert equation (|46[) by means of the same procedure as 
used for equation (|43|) . In addition, the integral on the left of equation (|46|) can also be calculated from the known spherically 
averaged protohalo density contrast profile. Thus, taking the ratios between the solutions for different j obtained from inversion 
of equation (|46[) . we can infer the eccentricity profiles for the seeds of typical haloes grown by PA. 

Once the eccentricity profiles, (e p ) p (r p ) and (e s ) p (r p ), are known, the typical halo shape and kinematics can be inferred 
following the steps described in Sections [3] and [4] 1) infer (through eqs. [25] and [38]) the eccentricity profiles e p (r) and e s (r) 
of the halo as a function of cr(r); 2) determine from them the mean squared density fluctuation profile for the halo (eq. [10] ^ 
as a function of <r(r); 3) infer the corresponding mean squared potential fluctuation profile (eqs [14] and [17]) as a function of 
<r(r); 4) solve the generalise Jeans equation (|42l) for the velocity dispersion profile er(r); and 5) determine (from eqs. [25] and 
[38]) the typical eccentricity profiles, e p (r) and e s (r), now using the explicit values of <r(r), and the anisotropy profile /3(r) 
(eq. [40] after deriving the mean squared potential fluctuation profile). 

Unfortunately, the accurate theoretical eccentricities so inferred could not be compared with the results of numerical 
simulations because there is so far no such typical profiles drawn from simulations. We just know some main trends of the 
typical shape of haloes. Thus, it is not worth at this stage to carry out such a complex derivation. In fact, the purpose of the 
present paper is not to infer accurate eccentricity and kinematic profiles for CDM haloes directly from the power-spectrum of 
the concordance model but rather to verify the validity of the model and try to understand the origin of the universal trends 
shown by such halo properties in numerical simulations. And these two objectives are much better obtained by means of the 
following simpler and comprehensive approximate procedure. 

Peaks are triaxial with a rather prolate shape (BBKS). As mentioned, their eccentricities depend on 5 p k and Rf through 
their curvature x in such a way that, the larger x, the more spherical is the peak. The quantity {{x)) increases progressively 
with increasing Rf over typical peak trajectories, meaning that peaks become increasingly spherical over those tracks, very 
slowly first and much more rapidly at the end. Thus, the same trend must be found in protohaloes for increasing r p as well 
as in typical haloes for increasing r. (It is true that the relation between the eccentricities in the seed and the halo also 
depends on cr(r), but this dependence can be neglected at small radii, while at large radii the decrease with increasing r p 
of protohalo eccentricities is so marked that the dependence on a(r) cannot reverse it.) Therefore, according to the present 
model, haloes should be rather prolate and approximately homologous at small and intermediate radii and ten d to become 
more spherical near the halo edge. This behaviour is fully consistent wit h the results of numerical simulations (e.g. Jing &: Sutol 
|2002| ; iBailin fc Steinmetjliooil ; 1 Allgood et alJbood : Istadel et afll200g| '). 



Moreover, we can have an idea on the overall values of halo eccentricities because, as mentioned, they are expected to be 
rather uniform except near the edge. The typical ratio between the major and minor axes (the major and intermediate axes) 
in peaks with low and moderate average curvature is ~ 1.7 (~ 1.3) (BBKS), implying typical values of the eccentricities e p 
(e s ) of ~ 0.81 (~ 0.64). Such typical values, independent of the filtering scale (they only show a moderate dependence on the 
curvature) were derived by BBKS for peaks in the old standard CDM model (i.e. in the Einstein-de Sitter cosmology), so they 
are valid for any CDM cosmology as all cosmologies converge for high enough redshifts as those corresponding to primordial 
peaks to the Einstein-de Sitter model. As U(r) is close to one with a small bias towards two at small r, the expected values of 
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Figure 3. Same as Figure 1 for asymptotic power-law regimes with the range of k and a indexes leading to only moderate values (not 
much greater than one) of Q/P as found in numerical simulations. 



e p (e s ) in haloes at small radii should be close to, perhaps a little larger than, those of peaks, say, ~ 0.9 (~ . 8). This is also 
consistent with the values found in numerical simulat i ons ( Frenk et al||l 988l; iBull ockl 1200 ^3; Ijing fc Sutc 20021; Springel et al, 
20041; iKasun fc Evrardll2005l ; iBailin fc Steinmet jlioOol ; I Allgood et al.lbood : Irlavashi et al.l l 



2007; 



Bctt et al 



23; E 



2007; Stadel et al 



2009). 



At large r, the dependence on a 2 (r) of the relation between the eccentricities in the halo and the seed becomes important. 
Yet, we can still determine the expected outer asymptotic behaviour of the eccentricities in the halo, taking only into account 
that, as mentioned, the object becomes more spherical at very large r. In that asymptotic regime, the power indexes of the 
mean squared density, squared potential and crossed density-potential fluctuation profiles coincide (see Sec. [3} and the more 
or less marked departure of isodensity and isopotential contours from spheres, given by the ratio Q/P, depends on the index 
k and the asymptotic logarithmic slope —a of the density profile (eq. [20]). For values of a approaching 3 as near the outer 
asymptotic regime of virialised haloes, the only possible negative values of n leading to outward-decreasing eccentricities with 
isopotential c ontours moderately more spherical than the isodensity contours {Q/P not much greater than one) as found in 
the literature (|Springel et alj|2004l ; IKasun fc Evrardll2005l ; lHavashi et alj|2007l ) are in the range —0.2 < k < (see the regions 
in dark blue, or quite strong grey if black and white, in Fig. [3}. Thus, the outer asymptotic logarithmic slope, n, is also 
severely constrained in this approximate treatment. 

And what about the typical halo kinematics? The typical mean squared potential fluctuation profile for CDM haloes 
in log-log should be approximately given by the natural spline between a uniform value of ~ 0.11, as implied by the above 
mentioned typical eccentricities, below about one hundredth the virial radius and a straight line with logarithmic slope k 
beyond about one tenth the virial radius. We will consider two values of n, —0.175 and —0.1, in order to better sample 
the allowed range —0.2 < n < and to see how robust the theoretical typical halo kinematics are against variations in the 
individual shapes of haloes. (Different triaxial shapes must give rise to different mean squared potential fluctuation profiles.) 
The two squared potential fluctuation profiles so built are shown in Figure [4] 

Once the form of the squared potential fluctuation profile has been fixed, we can calculate the ratio o\ (r)/a 2 (r) (eq. |40| ) 
and, from it, the anisotropy profile (eq. [41]). Figure [5] shows the velocity anisotropy profiles inferred from the two different 
squared potential fluctuation profiles. As can be seen, they are both in good agreement with the results of numerical simu- 
lations. Specifical l y, for n = —0.1, we find an anisotropy profile that closely follows the typical 'universal' law proposed by 
Hansen fc Stadel ( 20061 ) (the red line in the right panel). While, for n = —0.175, we find it closer to the anisotropy profile 



found bv lNavarro et al (|2010h in a realisation of a Milky Way mass halo (the red line in the left panel). The o nly noticeable 
disagreement in this latter case is that the theoretical profiles show no cutoff near the halo edge as found bv lNavarro et all 
But this is a very unusual feature (it does not follow the 'universal' trend) and likely reflects an incomplete relaxation of this 
particular simulated halo at those large radii. 

Using the previous theoretical velocity anisotropy profiles, we can solve the generalised Jeans equation (|42l) to find the 
corresponding velocity dispersion profiles and, using the typical spherically averaged halo density profile derived in SVMS, 
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Figure 4. Plausible typical mean squared potential fluctuation profiles according to our model. They connect in a smooth way an inner 
asymptotic uniform value consistent with those found in Figure [2] and outer asymptotic logarithmic slope k equal to —0.1 (solid line) 
and —0.175 (dashed line) within the allowed range according to Figure [51 
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Figure 5. Theoretical anisotropy profile derived from the rms potential fluctuation profiles shown in Figure[4](same lines). Left panel: as 
a func tion of radius to facilitate the comparison with the curve found for the Milky Way mass halo (red line) simulated bv lNavarro et all 
(|2010l : see their Fig. 10). Right pa nel: as a function of the logarithmic density slope to facilitate the comparison with the analytical 
formula by Hansen & Stadcl ( 2006) fitting the anisotropy profiles for simulated haloes (red line) . 
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Figure 6. Theoretical pscudo phase-space density profile resulting from the two anisotropy profiles shown in Figure [5] (same lines), 
compared to a straight line with slope —1.875 and arbitrary zero-point as found in simulations (black line). For comparison we also plot 
the theoretical profiles obtained for haloes with masses equal to 10 13 M Q (blue lines) and 10 14 M Q ( green lines). 

the corresponding pseudo phase-space density profiles, (p)(r)/a 3 (r). In Figure[6] we show the results obtained from the two 
mean squared potential fluctuation profiles above. In both cases, the theoretical pseudo phase-space density profile is very 
nearly a power-law with a logarithmic slope of —1.875, in full agreement with the results of numerical simulations. The two 
solutions essentially overlap; there is just a slight undulation at intermediate radii in the solution obtained from « = —0.175, 
which indicates that the spline carried on between the inner and outer asymptotes of the squared potential fluctuation profile 
is somewhat deficient in this case. Nonetheless, the fact that so distinct mean squared potential fluctuation profiles as adopted 
here give so similar solutions demonstrates that the pseudo phase-space density profile is very robust against variations in 
halo shape (see Sec. [6] for the possible origin of this behaviour). Note instead the clear dependence of the zero-point of the 
pseudo phase-space density profile on the mass of the halo predicted by the model. This cannot be compared with simulations 
as this particular dependence has not yet been analysed. It would be worth trying to confirm this prediction of the model. 

According to our results, the velocity anisotropy profile should show an important dependence on the shape of the system, 
which is consistent with the rather large scatter shown by this profile in numerical simulations. The velocity dispersion profile 
follows from the Jeans equation (|42l) for such an anisotropy profile and the density profile {p)(r). As the latter profile is 
independent of the shape of the object, it cannot balance in the Jeans equation the dependence of /3(r) on that property. 
Consequently, a(r) must also depend on the shape of the object. Yet, the pseudo phase-space density profile turns out to be 
very robust against variations in that property; there are just very slight undulations depending on it. Consequently, <j(r) 
must also be similarly robust. This is once again in agreement with the small scatter shown by this profile in simulated haloes. 

To end up we want to remark that we have made so far no reference to the CDM cosmology used; all the preceding 
(approximate) results hold, regardless of the cosmology provided only haloes grow hierarchically from the bottom-up (see 
SVMS). This is once again in agreement with the results of numerical simulations. Of course, the finer details of the predicted 
profiles will depend on the specific power-spectrum used, but this will not modify the main properties discussed here. 



6 DISCUSSION 

As just shown, the present model recovers the right form of the anisotropy and pseudo phase-space density profiles of simulated 
haloes. The origin of the former profile is clear: it arises from the tight relation (|40[) between the tangential velocity dispersion 
and the mean potential fluctuation profile and the typical form (essentially the well-constrained inner and outer asymptotic 
behaviours) of this latter profile. But the origin of the power-law-like pseudo phase-space density profile is less obvious. On 
the other hand, we have assumed so far that haloes grow by PA while they grow both through PA and major mergers. In this 
Section, we address these two important points. 
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6.1 Origin of the Pseudo Phase-Space Density Profile 

The form of the pseudo phase-space density profile of simulated haloes, recovered by the present model, coincides with that 
predicted by Bertschinger's (1985) model for collisionless spherically symmetric, self-similar, systems evolving by PA. This 
implies that the ultimate reason for such a coincidence cannot be neither the shape (spherically symmetric vs. triaxial) of the 
system nor the form (power-law vs. non-power-law) of the density profile. The only more fundamental aspect in the formation 
of virialised objects grown by PA directly or indirectly included in both models is the way the coarse-grained phase space 
density increases via phase mixing (entropy generation) when the system loses energy via shell-crossing and contracts. 

Bertschinger's model is dynamical and can properly follow the effects of shell-crossing, while the SVMS model is steady and 
focuses on the equilibrium state resulting from that process. However, the inside-out growth condition the SVMS model relies 
on directly follows from the lack of apocentre-crossing during virialisation, a condition that is also fulfilled by Bertschinger's 
model. Thus, despite the distinct approach, the effects of shell-crossing are similarly accounted for in the two models. 

Certainly, the power-law form of the pseudo phase-space density profile predicted by Bertschinger's model is tightly 
related to the self-similarity assumption, while the present model makes no such an assumption. However, PA is always close 
to self-similar. It is driven by gravitation, a scale-free force, and the initial conditions are also very approximately scale-free: 
the initial mass distribution is essentially uniform (with essentially the critical density) and the velocity field corresponds to 
the scale-free unperturbed Hubble-flow form. Thus, the solution of the present model at small z should not be very far from 
self-similar, except for the effects of the increasing departure from the Einstein-de Sitter universe owing to the specific density 
profile of the seed. In this sense, the predicted profiles should not be far from power-laws, particularly the pseudo phase-space 
density profile that appears to be so little sensitive to the details of the seed. Another fundamental assumption of Bertschinger's 
model is the radial infall of the system, while the SVMS model takes into account that the collapse and virialisation is actually 
non-radial. However, if the tangential velocities develop at the expense of the initial radial velocities without altering the total 
velocity dispersion, as assumed in equation (|40[> . the velocity dispersion that emerges as a consequence of shell-crossing should 
be the same in spherically symmetric as in triaxial systems. 

The situation is therefore as follows. The power-law form, with index equal to -1.875, of the pseudo phase-space density 
profile found by Bertschinger is the direct consequence of the shell-crossing produced in self-similar radial PA. Thus, the 
conditions that: i) PA is always quite close to self-similar and ii) the velocity dispersion generated by shell- crossing is not 
modified, according to the condition (|40|l . by the appearance of tangential velocities in non-radial PA, should guarantee that 
the evolution of the coarse-grained phase-space density found in Bertschinger's model is kept essentially unaltered in the 
general (non-striclty self-similar and non-radial) case of PA. This would explain why the present model recovers the pseudo 
phase-space density profile a la Bertschinger in agreement with the results of numerical simulations. Note that such a pseudo 
phase-space density profile does not necessarily implies that the density profile must also be roughly proportional to r -2,25 
as found by Bertschinger. This would be the case i f PA were st r ictly self-similar, so to end up with a power-law density 
profile, and, more importantly, radial. As stated bv iBertschingerl (1985) when comenting on his density profile p oc r -2 ' 25 : 



"Self-similar relaxation is not complete [as it pro ceeds along one dir ection only] and thus heuristically should not lead to a 
density profile as flat as p oc r~ 2 [as predicted bv iLvnden-Belll 1 1967 ) in the case of complete 3D relaxation]. It is uncertain 



to what extent the addition of angular momentum (i.e. non-radial orbits) will change this result,...". In other words, if PA is 
non-radial, than the power index of the density profile is expected to change, in agreement with the results of SVMS, despite 
the fact that the pseudo-phase space density profile, driven by the evolution of the phase-space density during virialisation in 
the absence of apocentre-crossing, should be kept roughly unaltered as a consequence of condition f|40[) . 

As mentioned, a key point in the present model, arising from the particular way shell-crossing proceeds in accreting 
virialised haloes, is their inside-out growth. Bertschinger's model is instead based on the self-similar evolution of such systems 
and does not care about the implications this has on the growth of the steady object. But, if our explanation of the power-law- 
like pseudo phase-space density profile is correct, the two models should essentially coincide (except for the symmetry and the 
strict self-similarity condition), meaning that Bertschinger's solution should also approximately satisfy the inside-out growth 
condition. This is indeed the case. As stated by Bertschinger, as a consequence of shell-crossing, particle orbits rapidly become 
approximately periodic, that is, the mass inside them stabilises and the inner system does not change anymore. In other words, 
the system develops a steady core and grows from the inside-out as new shells virialise. Of course, such a behaviour can only 
be approximately satisfied in Bertschinger's model. As a consequence of self- similarity, steadiness can only be strictly achieved 
in the limit of infinite time or vanishing radii; particle orbits can never become exactly periodic, they can only tend to become 
so (see Betschinger's Fig. 9 at the base of his comment regarding the inside-out growth). To obtain a fully steady solution, 
the self-similarity assumption must be replaced by the inside-out growth condition as in the present model. 



6.2 Major Mergers 

As has long bee n known, there is in radial PA a one-to-on e correspondence between the density profiles of virialised objects 
and their seeds ( Bertschinger 19851 : Del Popolo et al. 2000). The reason for this is that in radial PA, there is no phase- mixing 
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of particles turning around at the same time ( Bertschinger 19851 )) or, equivalently, there is no apocentre-crossing of particle 
orbits during virialisation. This result was extended in SVMS to non-radial PA. 

An important consequence of this result was that the spherically averaged density profile (p) (r) for a virialised halo does 
not allow one to tell whether or not it has suffered major merger^j. The reason for this is that, given a virialised halo with 
any arbitrary aggregation history, we can always think about one peak with appropriate density profile (or equivalently, the 
spherical energy distribution £ P (M)), leading by PA to a virialised object with identical spherically averaged density profile. 
Such a putative accreting seed of the halo really exists; it is the peak tracing the current evolution by smooth accretion of 
the halo according to the peak formalism. Thus, it is a normal peak contributing to the peak number density just as the halo 
associated with it is a normal halo contributing to the halo mass function. Therefore, the halo is indistinguishable from one 
evolved by PA, regardless of its real aggregation history. 

But what about the shape and kinematics of haloes? Do they allow one to tell between objects formed by PA and having 
suffered major mergers? The answer is always the same: major mergers go unnoticed in these properties as well. The proof is 
similar to that carried out in SVMS for the density profile. Using the one-to-one relations inferred in Section [3] thanks to the 
results of SMVS, given a halo with any arbitrary aggregation history, one can always think about one peak not only with the 
appropriate density profile, independent of the ellipsoidal shape of the peak, leading by PA to the spherically averaged profile 
of the halo (SVMS), but also with a shape also leading to its shape and kinematics. The shape of this putative accreting seed, 
coincides, by continuity, with that of the peak tracing the halo in the peak formalism (see SVMS) and is thus a normal peak 
contributing to the peak number density. As a consequence, the halo is indistinguishable from the one that would grow by 
PA from that seed not only regarding the spherically averaged density profile, but also the shape and kinematics. This result 
can be summarised by saying that major mergers go unnoticed in virialised objects regarding all those properties. 

As discussed in SVMS, the previous conclusion reflects the fact that virialisation is a real relaxation process yielding the 
memory loss of the halo past history. This obviously affects not only the inner structure but also the shape and kinematics of 
virialised objects. Certainly, the one-to-one correspondence between virialised haloes and their seeds i n PA implies that ther e 
is no full memory loss in PA due to virialisation, in agreement with the results of simulations (e.g. IVera-Ciro et alJl2011T ). 
However, the fact that it is impossible to reconstruct the seed of any simulated virialised halo by running the simulation 
backwards is evidence of the existence of a time arrow indicating the direction of virialisation. What makes both aspects 
consistent is the fact that, as just mentioned, the properties of virialised haloes do not allow one to tell whether they have 
suffered major mergers or they have evolved by PA. Therefore, we cannot unambiguously determine their initial conditions. 
In other words, we can reconstruct their putative accreting seeds, but these may or may not be their real (possibly multinode) 
seeds. 

An important corollary of this result is that the present model for the shape and kinematics of virialised objects formed 
by PA also holds for objects having suffered major mergers. One must simply infer the properties of the virialised object from 
those of its putative accreting seed, that is the peak tracing its current evolution by smooth accretion. 



7 SUMMARY AND CONCLUSIONS 

The model developed in SVMS for the inner structure (spherically averaged density profile) of virialised haloes has been 
extended to deal with their triaxial shape (eccentricity profiles) and kinematics (velocity dispersion and anisotropy profiles). 
To do this we have considered the simple scenario of haloes evolving by pure accretion. Under this assumption we have derived 
the shape and kinematics of the final objects from the shape of their seeds. The reason why dealing with the case of pure 
accretion is sufficient is that all virialised halos, even those having undergone major mergers, can be seen to arise in this 
way from the peak which, according to the peak formalism, traces their current evolution. The fact that in pure accretion 
the final objects keeps the memory of the initial conditions is a consequence of virialisation being achieved, in this case, 
through shell-crossing with no crossing of particle apocentres. This is not contradictory with the idea that virialisation is a 
real relaxation process. On the contrary, the fact that major mergers do not leave any particular imprint in the shape and 
kinematics of virialised objects implies that one cannot be sure about the real aggregation history or, equivalently, about the 
real seed of any given virialised object. 

Applied to virialised haloes in hierarchical cosmologies, the model allows one to make the link between the typical 
(mean) shape and kinematics of these objects and the power-spectrum of random Gaussian density perturbations. Following 
a simpler comprehensive approximate approach, we have shown that the overall shape and kinematics predicted by the model 
are consistent with the results of simulations of CDM cosmologies. In particular, the theoretical anisotropy profile has the 
same universal form as empirically found, with a substantial scatter due to the variety of peak (halo) shapes. Likewise, the 
theoretical pseudo phase-space density profile has very approximately a power-law form with the logarithmic slope of —1.875 
shown by simulated haloes. This profile is much more robust, that is, it does not essentially depend on the peak (halo) shape. 



We are of course presuming that the halo has completed relaxation after its last major merger. 
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The present model suggests that the origin of the power-law-like form of the pseudo phase-space density profile of virialised 
haloes formed by pure accretion is due to the particular evolution of the coarse-grained phase-space density during virialisation. 
This takes place through the phase mixing produced by shell-crossing in one direction with no crossing of particle apocentres. 
The fact that PA is always very nearly self-similar and that the velocity disp ersion is conserved w here a tangential velocity 
component develops in non-radial PA would explain why the solution found bv lBertschingerl (|1985T ) for self-similar, radial PA, 
with the same kind of phase mixing, is very nearly preserved. 
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APPENDIX A: SOLUTION METHOD FOR THE ECCENTRICITIES OF THE VIRIALISED OBJECT 
FROM THOSE OF THE PROTOHALO 



To solve equations (|25[) and (|38[) , this latter without the term in cr(r) 
and 



r ), or, equivalently (see eq. [10]), equations (|25ll 



(r) 



3[(l~ e 2 ) 2 (l- e 2 ) 2 + (l- e 2 ,) 2 + (l- e 2 ) 2 ] 
[(l- e 2 )(l-e 2 ) + (l-e 2 ) + (l-e 2 )]2 



(r P ), 



(Al) 



it is convenient to change the variables e p and e s into the quantities r c /r and r e /a. Given the definition of the eccentricities 
the relation between the old and new variables is 



(5) 6 = (l-4)(l-e 2 ) 



and 



(7)' 



27 



(l-egHl-e 2 ) 
[l + (l-e 2 ) + (l-e 2 )]3 



To invert them, we must first solve the biquadratic equation 

^HS>--S( 3 5-£)+ 2 =°. 

and then replace the solution into the equation 



(A2) 



(A3) 



(A4) 

Specifically, the only one solution of the biquadratic equation ()A3|) guaranteeing the positiveness of the eccentricities and the 



e P = x 1 - + 



el = 1 1 i- - -sHr I ~ 



1/2 



1/2 



In the new variables, equations ()25|) and (|A1|) take the form (see the relations [A2 



r e (r) _ r e (r ta ) 



and 



'■ta 



r rj(r) 
r e (r) a 4 (r) 



1 + £ U(r) 



5p_ 

(P) 



(r) 



1/2 r 



r b e (r) 



+ 3 



rl{r) 



r|(r) a 2 (r) 



(A5) 



(A6) 



(A7) 
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which can be readily solved for the quantities r e (r)/r and r/a(r). Then, the eccentricities e p (r) and e s (r) can be obtained by 
means of equations (|A5|) and ()A6[) . 
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